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Q\ • ABSTRACT 

in 

Context. About half of the baryons of the Universe are expected to be in the form of filaments of hot and low density intergalactic medium. Most 
of these baryons remain undetected even by the most advanced X-ray observatories which are limited in sensitivity to the diffuse low density 
medium. 

Aims. The Planck satellite has provided hundreds of detections of the hot gas in clusters of galaxies via the thermal Sunyaev-Zel'dovich (tSZ) 
. . . effect and is an ideal instrument for studying extended low density media through the tSZ effect. In this paper we use the Planck data to search for 
^> ■ signatures of a fraction of these missing baryons between pairs of galaxy clusters, 
r T ~j ' Methods. Cluster pairs are good candidates for searching for the hotter and denser phase of the intergalactic medium (which is more easily 
| observed through the SZ effect). Using an X-ray catalogue of clusters and the Planck data, we select physical pairs of clusters as candidates. Using 
, the Planck data we construct a local map of the tSZ effect centered on each pair of galaxy clusters. ROSAT data is used to construct X-ray maps of 
these pairs. After having modelled and subtracted the tSZ effect and X-ray emission for each cluster in the pair we study the residuals on both the 
SZ and X-ray maps. 

Results. For the merging cluster pair A399-A401 we observe a significant tSZ effect signal in the intercluster region beyond the virial radii of 
the clusters. A joint X-ray SZ analysis allows us to constrain the temperature and density of this intercluster medium. We obtain a temperature of 
kT = 7.1 ± 0.9 keV (consistent with previous estimates) and a baryon density of (3.7 ± 0.2) x 10~ 4 cm 3 . 
Conclusions. The Planck satellite mission has provided the first SZ detection of the hot and diffuse intercluster gas. 

Key words. galaxies:clusters:general 



1. Introduction 

A sizeable fraction of the baryons of the Universe are expected to 

be in the form of the WHIM (Warm-Hot intergalactic medium) 

*** Corresponding author: G. Hurier hurier@lpsc.in2p3.fr and remain undetected at low redshifts. The WHIM is expected 
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to exist mostly in filaments but also around and between mas- 
sive clusters. These missing baryons are supposed to be in a 
low-density, low-temperature phase (over-densities between 5 
to 200 times the crit ical density and T = 10 5 - 10 7 K, 
ICen&Ostrikeri (fl999)) making the amount of X-rays produced 
by the WHIM too small to be detected with current X-ray fa- 
cilities. By contrast, their detection cou ld be possible via the 
Sunyaev-Zel'dovich (SZ hereafter) effect dSunvaev & Zeldovichl 
1 19721) produced by the inverse Compton scattering between the 
Cosmic Microwave Background (CMB) photons and the elec- 
trons of the WHIM. As the SZ effect is proportional to the elec- 
tron pressure in the medium, low density and low temperature 
regions can be detected provided their integrated signal is large 
enough. Planck's relatively poor resolution becomes an advan- 
tage in this situation since it permits scanning of wide regions of 
the sky that can later be integrated to increase the signal-to-noise 
ratio of the diffuse (but intrinsically large-scale) SZ signal. 

The full-sky coverage and wide frequency range of the 
Planck satellite mission makes it possible to produce reli- 
able ma ps of the tSZ emission dPlanck Coll aboratio n et al.l 
1201 lbldlla) . In particular Planck is better suited than ground ex- 
periments to detecting diffuse SZ signals, like the WHIM, which 
can extend over relatively large angular scales. Ground exper- 
iments can be affected at large angular scales by atmospheric 
fluctuations that need to be removed. This removal process can 
distort the modes that include the large angular scale signals. 
Planck data does not suffer from these limitations and can use 
its relatively poor angular resolution (when compared to some 
ground experiments) to its advantage. Indeed, diffuse low sur- 
face brightness objects can be resolved and detected by Planck. 
Finally, the wide frequency coverage and extremely high sensi- 
tivity of Planck allows for detailed foreground (and CMB) re- 
moval that otherwise would overwhelm the weak signal of the 
WHIM. 

The gas around clusters is expected to be hotter and denser 
than the WHIM in filaments, making direct detection of the 
cluster gas more likely. In addition, the increase of pressure 
caused by the merging process enhances the SZ signal, making 
it easier to detect the gas between pairs of interacting clusters. 
In the process of hierarchical formation clusters assemble via 
continuous accretion and merger events. Therefore, the bridge 
of intercluster matter between them is expected to be of higher 
density, temperature and thus thermal pres sure than the aver - 
age WHIM matter found in cosmic filaments dDolag et al. 2006). 

The Planck satellite dPlanck Collaboration et alJ d2011alV ) 
has the potential to detect these filamentary structures directly 
via the SZ effect. Suitable targets for Planck are close objects 
which subtend large solid angles and therefore have large 
integrated SZ fluxes. Alternatively, regions between mergers 
(filaments between pairs of clusters) or extremely deep grav- 
itational wells (superclusters like the Shapley) will contain 
diffuse gas with increased pressure which could be detected 
by Planck. For this work, we concentrate on searching for 
diffuse filamentary-like structure between pairs of merging 
clusters. We used th e MCXC catalogue of clusters of galaxies 
(Piffa retti et al.ll201ll) and the Planck data to select a sample of 
pairs of merging clusters in order to study the properties of the 
gas in the intercluster region. 

Indirect WHIM detections have been claimed through ab - 
sorption lines in the X-ray (and UV) band (Ric hter et al] l2008). 
There is also evidence of filamentary structure in the inter- 
cluster region from X-ray observations of several well known 



merging clus ter pairs like A222- A223 dWerner et alJ |2008), 
A399-A401 dSakelliou & Ponmanl 12004 and A3391-A3395 
(Tittlev & Henriksen 2001), and from the double cluster A1758 
(burretetal]l201lh . The pairs of clusters A3391-A3395 (sep- 
arated by about 50' on the sky and at red shifts z=0.051 and 
z=0.057 respectively dTittlev & Henriksenl20 0D) and more spe- 
cially, A399-A401 (separated by about 40' on the sky and at 
redshifts z=0.0724 and z=0.0737 respectively) are of particular 
interest for the purpose of this paper, given their geometry and 
angular separation. This is sufficient to allow Planck to resolve 
the individual cluster components. 

For A399-A401, earlier observations show an excess of 
X-ray emission above th e background level in the in tercluster 
region. Using XMM data, Sakelli ou & Ponmanl d2004l) obtained 
best fitting models in the intercluster region which indicated 
such an excess. Both clusters are classified as non-co ol-core 
clusters and show weak radio halos dMurgia et alj|2010l) . These 
two facts could be an indication o f a past interaction between 
the two clusters. iFuiita et alJ dl996l) analyzed ASCA data of the 
intercluster region and found a relatively high temperature in 
this region. They suggest ed a pre-merger sc enario but did not 
rule out a past interaction. iFabian et all d 19971) used HRJ ROSAT 
data and found a prominent linear feature in A399 pointing 
towards A401. They suggested that this coul d be evidence of a 
past interaction. Using Suzaku observations, dFuiita et al.l f2008) 
found that the intercluster region has a relatively high metallicity 
of 0.2 solar. These works estimate that the filamentary bridge 
has an electron density of n P ~ 10" 4 cm" 3 dFuiita et al.l 1 1996b 
ISakelliou & Ponmanll2004lFuiita et alj|2008t) . 

In this paper we concentrate on pairs of merging clusters 
including A399-A401 and we study the physical properties of 
the gas in the intercluster region via a combined analysis of the 
tSZ effect and the X-ray emission. The paper is organized as 
follows. Section [2] gives a brief description of the Planck data 
used for this study. In Section [3] we describe the selection pro- 
cedure used to identify the most suitable pairs of clusters for the 
analysis. We search for pairs of clusters for which the contribu- 
tion of the SZ effect to the signal is significant in the intercluster 
medium. Section|4]describes the X-ray ROSAT observations for 
the selected pairs of clusters. In section|5]we model the SZ and 
X-ray emission from the clusters assuming spherical symmetry 
and subtract them from the data. In section [6] the SZ and X-ray 
residuals are fitted to a simplified filament model to characterize 
the physical properties of the intercluster region. Section [8] dis- 
cusses our main results focusing on the limitations imposed by 
the cluster spherical symmetry assumption and alternative non- 
symmetric scenarios. Finally, we conclude in Section[9] 

2. Planck data 

Planck dTauber et al.ll2010t iPlanck Collaboration 1 1201 ih is the 
third generation space mission to measure the anisotropy of the 
cosmic microwave background (CMB). It observes the sky in 
nine frequency bands covering 30-857 GHz with high sensitiv- 
ity and angular r esolution from 31' to 5'. The Low Frequency 
Instrument LFI ; (Mand olesi et alj|2010 ; Bersa nelli et alJ 12010; 
iMennella et al.ll201 lh covers the 30. 44, and 70 GHz bands with 
amplifiers cooled to 20 K. The High Frequency Instru ment (HFI; 
lLamarre et al.l 1201 Ot IPlanck HFI C ore TearrJ |201 lal) covers the 
100, 143, 217, 353, 545, and 857 GHz bands with bolome- 
ters cooled to 0. 1 K. Polarization is measured in all bu t the 
highest two bands dLeahv et alJ 1201 Ot iRosset etal] 120101) . A 
combination of radiative cooling and three mechanical coolers 
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Table 1. Main physical parameters of the selected pairs of clusters. Values quoted with NA are unknown. 
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produces the temperatures nee ded for the detectors and optics 
dPlanck Collaboration" Two Data Processing Centers 

DPCs) check and calibrate the data and make maps of the sky 
Planck HFI Core Team! 1201 lbh IZacchei et all 1201 ll) . Planck's 
sensitivity, angular resolution, and frequency coverage make it a 
powerful instrument for galactic and extragalactic astrophysics 
as well as cosmology. Early astrophysics results are given in 
Planck Collaboration VIII-XXVI 201 1, based on data taken be- 
tween 13 August 2009 and 7 June 2010. Intermediate astro- 
physics results are now being presented in a series of papers 
based on data taken between 13 August 2009 and 27 November 
2010. 

This paper is based on Planck's first 14 month survey 
mission. Th e whole sky has been covered more t han two times. 
We refer to iPlanck HFI Core Team etaP (l201lh : IZacchei et all 
(1201 11) for the generic scheme of TOI processing and map- 
making, as well as for the technical characteristics of the maps 
used. This work is based on the nominal survey full-sky maps 
in the nine Planc k frequency bands provided in HEALPIX 
dGorski et al.ll2005l) with nside=2048 and full resolution. An er- 
ror map is associated with each frequency band and is obtained 
from the difference of the first half and second half of the Planck 
rings for a given position of the satellite. The resulting maps 
are basically free from astrophysical emission, but they are a 
good representation of the statistical instrumental noise and 
systematic error. We adopted circular Gaussian beam patterns 
with FWHM values of 32.6, 27.0, 13.0, 9.88, 7.18, 4.87, 4.65, 
4.72 and 4.39 'for channel frequencies 30, 44, 70, 100, 143, 
217, 353, 545 and 857 GHz, respectively. The uncertainties in 
flux measurements due to beam corrections, map calibrations 
and uncertainties in ban dpasses are expected to be smal l, as 
discussed extensively in (Planck Collaboration et al. 201 lb d.c) 



3. Pairs of merging clusters in Planck 

In the region trapped between interacting pairs of clusters we 
expect a hotter and denser phase of the WHIM (see Section []]) 
that might produce enough SZ emission to be detected by the 
Planck satellite. The combined full sky capabilities, lack of 
atmospheric fluctuations, and wide frequency coverage makes 
Planck a unique instrument to study these peculiar objects. 

To identify pairs of merger clusters in the Planck maps 
we use the MCXC catalogue of clusters Piffar etti et alJ (1201 lh 
that provides the main physical parameters of X-ray detected 
clusters of galaxies (position, redshift, angular extension 6*500 
and mass Msoo). We select a sample of cluster pairs for which 
their difference in redshift is smaller than 0.01 and their angular 
distance is between 10 and 120'. For all the selected pairs we 
construct maps of the tSZ (tSZ hereafter) emission from the 
Planck HFI frequency maps at a resolution of 7.1 8' using dif- 
feren t compon ent separation techniqu es: MILCA (H urier et al. 
2010), NILC dRemazeilles et alJl201lh and GMCA (Bobi n et al. 
2008). A detailed discussion of the relative performance of these 
compo nent separation tech niques can be found in a companion 
paper (Collaborationl l2012h . In what follows we will take as a 
reference the MILCA maps to simplify the discussion. The tSZ 
maps are centered on the barycentre of the X-ray position of the 
two clusters and extend up to five times the distance between the 
two clusters (as given in the MCXC catalogue). The error dis- 
tributions on the final SZ maps ar e computed using the Planck 
Jack-k nife maps described in iPlanck HFI Core Team et al] 
(l20Tlh . 

Using these tSZ maps, we finally select those pairs of 
clusters for which at least one of the clusters has a signal- 
to-noise ratio larger than five. Table Q] shows the position (in 
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a) tSZ for A0399— A040! 



b) tSZ for A202S— A203 
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Fig. 1. MILCA maps of the Compton parameter y x 10 6 for the selected pairs of clusters. From left to right and from top to bottom 
we show the pairs of clusters a) A0399-A0401, b) A2029-A2033, c) A2147-A2152, d) A2256-A2271, e) MKW 3s-A2063, f) 
A3391-A3395 and g) A0209-A0222. 



Galactic coordinates), the redshift, and the angular diameter (in 
05oo units) of the two clusters as well as the angular distance 
between them (in arcminutes) for each of the selected pairs 
of clusters. We have divided the selected clusters into three 
different sets (separated by horizontal lines in Table [TJ. The 
first set corresponds to the pairs of clusters for which both 
clusters are significantly detected in SZ by Planck. The second 
set corresponds to the pairs for which one of the clusters 
is only marginally detected in SZ by Planck. The third set 
consists of pairs of clusters separated by less than 30 '. After 
a first assessment we concludee that for the last two sets the 
modelling of the clusters in the pair is too complex, given the 
lower significance of the data. This makes it difficult to extract 
reliable information on the intercluster region, therefore we do 
not consider them further in the analysis presented in the next 
sections. 



Figure Q] shows maps of the Compton parameter y x 10 6 for 
the six pairs of clusters in the first group in Table Q] (in de- 
creasing order of y from a) to f) ). These maps are rotated so 
that the virtual line connecting the centre of the clusters is hor- 
izontal and at the middle of the figures. The direction defined 
by this line will be called longitudinal hereafter, and the direc- 
tion perpendicular to it will be called radial. The pixel size in 
the longitudinal and radial directions is defined by the angular 
distance between the centre of the two clusters, d, as follows: 
Qpix = 5 x d/60. From visual inspection of these maps we infer 
that there might be significant SZ emission in the intercluster re- 
gion of the pairs of clusters A399-A401, A2147-A2152, MKW 
3s-A2063 and A3391-A3395. To confirm these findings we have 
computed 1 -D longitudinal profiles by stacking the above maps 
in the radial direction. Figure [2] shows the 1-D profiles for the 
six pairs of clusters (notice that the errors are correlated). The 
figure also shows the residuals after subtracting an estimate of 
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the tSZ emission from the clusters. In the intercluster region we 
simply use a symmetric interpolation (in red in the figure) of 
the external profile of each of the clusters. When there was an 
obvious contribution from extra clusters (see for example the 
shoulders on the A3391-A3395 pair (Figure 2e) and the peak 
on the A2029-2033 pair (Figure 2b)) we excluded the region 
affected from the analysis. The error bars on the residuals are 
increased by a factor \j2 in the outskirts and V3 in the inter- 
cluster region to account for the interpolation, symmetrization 
and subtraction. We observe clear extra emission for two of the 
pairs: A399-A401 (Figure 2a) and A3391-A3395 (Figure 2e). In 
the intercluster region, x 2 f° r the nu ll hypothesis (no extra tSZ 
signal) is 20 for A399-A401 and 174 for A3391-A3395, for 8 
degrees of freedom (11 data samples minus 3 parameters). We 
have also found large x 1 values for the A2029-A2033 (Figure 
2b) and A2256-A2271 (Figure 2d) pairs of clusters. However, 
for these two pairs we observe a deficit of tSZ signal in the inter- 
cluster region that is induced either by an over-estimation of the 
tSZ emission of the clusters themselves or by their asymmetric 
geometry. Thus, we only consider the A399-A401 and A3391- 
A3395 pairs of clusters for further analysis as they show the most 
prominent tSZ signal in the intercluster region. To check that the 
excess of tSZ effect observed can not be explained by contami- 
nation by foreground emissions (Galactic) we have estimated the 
dust temperature using Planck and IRAS data and a simple mod- 
ified blackbody SED that we fit to the Planck+IRAS data. No 
large deviations in the temperature that could compromise the 
component separation process are observed in the field of view. 

4. X-ray data for the selected pairs 

To improve the quality of our analysis we complement the 
Planck data with X-ray observations retrieved from the ROSAT 
archive. Our choice of using ROSAT/PSPC is motivated by 
i) its ability to detect the faint surfac e brightness emission 
(i.e the cluster em ission at large radii dVikhlinin et alj [T"999t 
lEckert et alJ 201 1[) and t he intercluster region in superclusters 
(iKull & Bohringeri ll999)). ii) its large field of view (~2 deg 2 ) 
and iii) the low instrumental background. These factors make 
ROSAT/PSPC the X-ray instrument which is most sensitive to 
diffuse low surface brightness emission. On the negative side, 
its limited bandpass and poor spectral resolution is not optimal 
for measuring plasma temperatures. 

We use the ROSAT Exten ded Source Analysis Software 



(ESAS, ISnowden et all (11994))) for data reduction, following 
the procedure in lEckert et al.l d201 ll) . We produce images of 
the counts in the R37 ROSAT band, i.e. in the energy range 
0.42 - 2.01 keV, and we correct for vignetting effects through 
the exposure map p roduced with the ESAS software. Following 
lEckert et alj d201 lb . we detect point sources up to a constant 
flux threshold, in order to resolve the same fraction of the 
CXB on the entire field of view ( FOV). We estimate the 
sky background components in the fraction of the FOV not 
contaminated by cluster emission (r > 1.3r2oo) and then subtract 
this constant value from the image. Error bars for the count rate 
were estimated by propagating Poissonian errors of the original 
source and background counts images. 

Regarding the ROSAT/PSPC point spread function (PSF), 
this strongly depends on the position in the FOV, ranging from 
15" on axis to 2' in the outer parts of the FOV. We consider 
an average PSF representative of our data sets as computed 
assuming a position 1' off axis and having 1 keV energy. 



We have checked that variations within reasonable limits 
(E = 0.5 - 2keV and OffAxis = 0-10 arcmin) do not introduce 
significant changes in our results. In order to reduce statistical 
noise we smooth the X-ray data with a 2' Gaussian kernel. The 
scale of this smoothing is significantly smaller than Planck's 
angular resolution so it does not compromise the results derived 
from X-ray data. This also helps to make the uncertainties in the 
ROSAT PSF a second order effect. 

Since we use PSPC archive data and we cover regions that 
might be larger than the FOV of ROSAT, we need to combine 
different pointings into a single mosaic. The pair A399-A401 
(Figure 3a) is contained in a single PSPC pointing so in this case 
it is not in principle necessary to perform a mosaic. However, 
we still combine the neighboring pointing rp800182n00 and 
rp80235n00, to increase the statistics. The total combined ex- 
posure time is ~ 14 ks. In the case of A3391-A3395 (Figure 
3b), the larger separation in the sky between the clusters forces 
us to combine at least two pointings. We combined the two avail- 
able observations for this pair: rp800079n00 (2.5 ks centered on 
A3395) and rp800080n00 (6 ks centered on A3391). The final 
X-ray maps in counts per second are presented in FigurefJ] 

5. Modelling of the tSZ and X-ray emission from the 
clusters 

One of the key issues of the analysis presented in this paper 
is the estimation of the tSZ and X-ray emission of the clusters 
themselves since this estimate is later used to remove the contri- 
bution of the clusters from the intercluster region. In this section 
we present a detailed description of the modelling of the clus- 
ter emission for the A399-A401 and A3391-A3395 cluster pairs. 



5. 1 . tSZ and X-ray emission from clusters 

When CMB photons cross a galaxy cluster there is the probabil- 
ity that some of the photons will interact with the free electrons 
in the hot plasma through inverse Compton scattering. The tem- 
perature decrement (increment for frequencies v > 217 GHz) 
observed in the direction 6 can be described as; 



AT(ff) 



= T„ j n , 



(l)T(l)dl 



(1) 



where T contains all the relevant constants including the fre- 
quency dependence (g x = x(e x + l)/(e v - 1) -4 with x = hv/kT), 
n e is the electron density and T is the electron temperature. 
The integral is performed along the line of sight. The quantity 
P([) - n e (l)T([) is often referred as the pressure produced by the 
plasma of thermal electrons along the line of sight. 

We can also consider the thermal X-ray emission from the 
hot ionized plasma including bremsstrahlung and line and re- 
combination emission from metals given by 



SM =S 



J n 2 e A(T)dl 



Di(z) 2 



(2) 



where Di(z) is the luminosity distance. The quantity S„ contains 
all the relevant constants and corrections (including the band 
and k-corrections). Combining X-ray and SZE observations it 
is in principle possible to break the degeneracy between differ- 
ent models due to their different dependency on T and specially 
on n e . 



5 







Fig. 2. Like in figure Q] From left to right and from top to bottom, 1-D tSZ longitudinal profiles and residuals after subtracting the 
contribution from the clusters (see text for details) for the pairs of clusters a) A399-A401, b) A2029-A2033, c) A2147-A2152, d) 
A2256-A2271, e) A3391-A3395, f) MKW 3s-A2063, and g) A0209-A222. 
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Fig.3. Maps of the X-ray emission (in counts per second) for the cluster pairs a) A3391-A3395 and b)A399-A401. 



5.2. Pressure profile models 

The tSZ and X-ray emission of clusters can be computed from 
the density and temperature profiles of the thermal electrons in 
the cluster. Below we describe the models of the electron pres- 
sure profiles for clusters that will be used to estimate and subtract 
the cluster contributions from the pairs of clusters. 



p-model 

The isothermal yS-model is historically the most widely used in 
the literature in the context of galaxy clusters. Even if cluster 
observations have shown that this model is over-simplified and 
unable to describe the details of the inter cluster medium (ICM) 
distribution, it allows us to use a limited number of parameters 
to produce a useful first approximation of the radial behaviour 
of the gas. In fact, we are not interested in the details of the 
single objects but in removing the main signal coming from the 
two clusters in the system in order to study the SZ excess signal 
observed between them by Planck. In the isothermal /?-model, 
the electron temperature is assumed to be constant and its density 
is given by; 



n e (r) = n e0 1 + 



(3) 



where n e Q is the central electron density and r c is the core radius. 



Generalized Navarro Frenk and White (GNFW) pressure pro- 
file 

Given the characteristic pressure Pgoo, defined as 



P 500 = 1.65 • HT^(z) 



M 5O0 



3 ■ lO 14 /!^ 



Mr 



2/3 



h^keVcnf, (4) 



the di mensionless universal (scaled) pressure profile can be writ- 
ten as ( Nae ai et alJl2007t lArnaud et alJl2010l) 



P(x) = 



Vo 



in which x - r/R^oo, c^qo = Psoo/rs, and y, a and (3 are re- 
spectively the central (r « r s ), intermediate (r ~ r s ) and outer 
(r » r„) slopes. The radial pressure can be the expressed as 



P(r) = P 500 F(x), 



(6) 



with [ Po , csoo, y,o-,ffl = \SA03h^ n , 1.117, 0.3081, 1 .0510, 5.4905], 
as derived by lArnaud et al.l (l2010t) . For the sake of simplicity, 
we redefine the normalization as Po = P500 * Po, also in units of 
keVcm 3 . 

As a result of the anti-correlation of the density and tem- 
perature profiles, at cluster cores the pressure exhibits less scat- 
ter than the electron density and temperature, and so is better 
suited to the definition of a universal profile. However, in the 
core region, t he disp ersion about the average profile found by 
lArnaud et all d2010li is still significant (~ 80 % at 0.03 P 50 o), 
while it becomes less than 30% beyond 0.2 Psoo- 



5.3. Fitting the tSZ and X-ray emissions from the clusters 

We use the /3 and GNFW models to fit the tSZ and X-ray 
emissions of the clusters in the A399-A401 and A3391-A3395 
pairs. In the case of the /3-model we vary the central density, 
n , the slope /3 and the core radius r c of each cluster. For the 
GNFW model we have considered two different sets of models: 
GNFW1 and GNFW2. For the first one, we fix j3 = 3.5 and fit 
the ot her parameters. For t he second one, C500, J and a are fixed 
to the lArnaud et al.1 d2010h values, leaving Po and /? free to vary. 

We fit each cluster individually excluding the intercluster 
region from the fit. For the tSZ emission we integrate the 
pressure profile along the line of sight up to 5 x P500 (unless 
stated otherwise) following equation Q] The resulting map 
is convolved with a Gaussian beam of 7.18' to match the 
resolution of the MILCA SZ map. For the X-ray emission 
we consider an isothermal scenario and compute the electron 
density from the pressure profile. We then use equation [2] 
to compute the surface brightness of the clusters. An X-ray 
map in counts per s econd is ob t ained from t he latter using 



(c 5 oo*) r [1 + (csooxT] 



a-\(J3-y)/a 



map in counts per s econd is OD t ainea rrom t ne l atter using 
( } the MEKAL model dMewe et al.l 1 19851 119861: iKaastral 119921: 



7 



Planck Collaboration: P. A. R. Ade et al.: SZ in merging clusters 



4iA 





n 




Fig. 4. a) Planck tSZ Compton parameter map; b) residuals after subtraction of the GNFW2 model; c) residuals after subtraction of 
the GNFW1 model; and d) residuals after subtraction of the /3 model of the clusters as described in the text. 
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Fig. 5. X-ray data and residuals of the best-fit cluster model for the A399 and A401 pair of clusters, a) ROSAT map; b) residuals for 
the GNFW2 model; c) residuals for the GNFW1 model; and d) residuals for the (3 model. 



iedahletal 1I1995I) for cl uster emission and the WABS model 
Morrison & McCammo^ 119831) for absorption of the neutral 
hydrogen along the line-of-sight. The MEKAL model is a 
function of the square of the electron density, the temperature 
and the redshi ft of the cluster. The WABS model uses colum n 
density maps ( Dick ev& LockmarJ 1 1 990t iKalberla et al. 2005). 
We convolve the final X-ray map with a Gaussian beam of 2' to 
match the resolution of the 2' resolution (degraded) ROSAT map. 

We have performed different fits including X-ray data only, 
SZ data only and both X-ray and SZ data. The fits are performed 
directly on the 2-D maps presented above and correlated errors 
are accounted for in the likelihood analysis. For the joint tSZ 
and X-ray fits the likelihood functions are normalized by their 



volume and then multiplied to obtain the best-fit parameters and 
errors. 



Results on the A399-A401 pair 

We converted the density profile into ROSAT PSPC cou nt rates, 
using an absorbed MEKAL model within XSPEC dArnaudl 
1996) . The conversion r a te is c omputed using the temperatures 
in ISakelliou & Ponmanl d2004l) . a metal abundance Z = 0.2 
solar, the Galactic colu mn density in the directi on of the pair 
nH = 1.09 x 10 21 cm" 2 dDickev &Lockmanll 19901) and the red- 
shifts z = 0.0737 for A401 and z = 0.0724 for A399. We fixed 
the truncation radius in the 3-D integrals to the r^oo value in 
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Fig. 6. 1-D a) tSZ and b) X-ray profiles and residuals for the best- fit models of the A399 and A401 clusters. The red curve corre- 
sponds to the GNFW2 model. The dotted lines are the model from Sakelli ou & Ponmanl d2004l) . 

Table 2. Best-fit parameters for the pressure profile /3-model. Notice that for this analysis KT and R max were fixed to the values 
shown. 



Cluster 


P 


r c (kpc) 


n (cm 3) 


KT (keV) 


Rmax (kpc) 


A399 


0.928 


432.8 


2.37 x 10" 3 


7.23 


2160 


A401 


0.928 


298.4 


5.50 x 10~ 3 


8.70 


2340 


A3391 


0.620 


213.0 


2.61 x 10~ 3 


6.00 


1500 


A3395SW 


0.673 


328.0 


1.17 x 10~ 3 


4.80 


1400 


A3395E 


0.726 


356.6 


1.17 x 10" 3 


5.00 


1200 



Table 3. Best-fit parameters for the GNFW1 pressure profile model. Notice that for this analysis j3 and KT were fixed to the values 
shown. 



Cluster 


C500 


7 


a 


P 


r 5 oo (Mpc) 


P (keV / cm 3 ) 


kT (keV) 


A399 


1.80 


0.000 


1.25 


3.5 


1.11 


2.04 xl0~ 2 


7.23 


A401 


2.58 


0.016 


1.59 


3.5 


1.24 


3.76 xlO 2 


8.47 


A3391 


0.69 


0.033 


0.72 


3.5 


0.90 


3.08 xlO- 2 


6.00 


A3395SW 


2.65 


0.066 


2.00 


3.5 


1.10 


0.70 xlO" 2 


4.80 


A3395E 


0.3 


0.000 


0.58 


3.5 


0.90 


1.89 xlO 2 


5.00 



Sakelliou & Ponman (2004) (this choice is not critical provided 
the radius is reasonably large). 

Before model fitting, we use three different methods to 
estimate the background or zero level of the tSZ Compton 
parameter maps. The results are consistent between the three 
methods. The first method computes the mean signal in an area 
around the cluster pairs that excludes the pairs. This region 
does not show any other source that needs to be masked out. 
The second method computes the mean signal of the averaged 
profiles of the clusters beyond R200 and in the directions 
opposite to the other cluster. The third method computes a 
histogram of the background region defined in the first method 
and finds the median of the distribution. We adopt method one 
for the final results. 

Tables|2][3]and|4]present the best-fit parameters of the A399 
and A401 clusters for the fj, GNFW1 and GNFW2 pressure 
profile models. The residual tSZ maps for these best-fit cluster 
models are shown in Figure |4] We notice that for the GNFW2 
pressure profile model (Figure 4d) the derived temperatures are 
lower than those derived from X-ray only data and in particular 
with t he results from XMM-Newton by ISakelliou & Ponmanl 
(2004j ). For the other two models the temperatures were fixed 
to the ISakelliou & Ponmanl d2004j) value. We show the Planck 



tSZ map (Figure 4a), and the residuals after subtraction of 
the best-fit f3 model (Figure 4b) and of the best-fit GNFW 
model for the two sets described above (Figure 4c and 4d). 
In the residuals we clearly observe an excess of tSZ emission 
with respect to the background. This can also be observed in 
Figure |6^) where we present the 1-D tSZ longitudinal profile 
an d residuals for the GNFW2 model described before and for 
the ISakelliou & Ponmanl (12004) best-fit model . 

In Figure [5] we show the ROSAT X-ray maps for the A399 
and A401 clusters and the residuals after subtracting the /3, 
GNFW1 and GNFW2 models. We observe that the models 
slightly under-estimate the signal at the centre of the clus- 
ters. This is also seen in Figure |6j)) where we present the 1- 
D x-ray longitudinal profile and residuals for the three pres- 
sure profile models. A si milar behavior was observed when fit- 
ting XMM-Newton data inlSakelliou & Ponmanl d2004l) . Quoting 
Sakelliou & Ponman (2004) neither of the two central galaxies 
is known to host an active nucleus, whose presence could be in- 
voked to explain the requirement for an extra, but small, central 
component. For illustration we als o plot (dotted lines) the best - 
fit model of the clusters found by ISakelliou & Ponmanl (12004). 
The latter model clearly over-predicts the tSZ cluster emission. 
However this discrepancy is not surprising, given the different 
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Table 4. Best-fit parameters for the GNFW2 pressure profile model. Notice that for this model C500, y, a and rsoo were fixed to the 
values shown. 



Cluster 


c 500 


r 





P 


r 50 o (Mpc) 


P (keV cm 3 ) 


kT (keV) 


A399 


1.18 


0.308 


1.05 


3.5 


1.11 


0.97 xl0~ 2 


5.0 


A401 


1.18 


0.308 


1.05 


5.0 


1.24 


2.40 xl0~ 2 


6.5 


A3391 


1.18 


0.308 


1.05 


3.5 


0.89 


0.75 xl0~ 2 


6.0 


A3395SW 


1.18 


0.308 


1.05 


3.5 


0.93 


0.40 xlO- 2 


4.8 


A3395E 


1.18 


0.308 


1.05 


3.5 


0.93 


0.40 xlO 2 


5.0 



dependencies of the SZ and X-ray signa ls on the ICM densit y 
and temperature. In fact, as pointed out by Hallm arTet al.l(l2007l) . 
the inadequacy of the isothermal assumptions affects the param- 
eter determination performed on SZ and X-ray data in a different 
way, therefore X and SZ profiles cannot be represented prop- 
erly using the same parameter values for an isothermal yS-model. 
The different sensitivity to the details of temperature and gas 
density distribution is still relevant even when comparing inte- 
grated quantities such as Y, the integrated Comptonizat ion pa - 
rameter, and the proxy Yx - kBTM gas , dKravtsov et al J 12006). 
whose ratio ha s been found to be smaller than 1 by several dif- 
ferent authors ( Arnaud et al. (2010); Planck Collaboration et al. 
(201 lp); lAndersson & SPT Collaboration! (120101) ; iRozo et al.1 
(2012)). This is in agreement with our findings and confirms the 
importance of a joint X-ray/SZ analysis, able to break the degen- 
eracy between models. Nevertheless, despite the uncertainties in 
the central part of the clusters, we observe an excess of X-ray and 
tSZ signal in the intercluster region. This excess is discussed in 
section[6] 



plex model is necessary to subtract the cluster components. We 
can also notice in the figure the contribution from a radio point- 
source that appears as a decrement in the reconstructed Compton 
parameter map. The area contaminated by this point source was 
excluded from the analysis. Due to the uncertainties in the model 
and in the foreground subtraction (radio point source contribu- 
tion), the analysis of the intercluster region gives unreliable con- 
straints on the density and temperature of this region and it is not 
considered in the following analysis. 



6. Analysis of the intercluster residuals 

6.1. Models 

To model the intercluster region of the pair of clusters we con- 
sider either an extra background cluster or a filament-like struc- 
ture described by a parametric model. 



Results on the A3391-A3395 pair 

The modelling of the A3391 and A3395 pair is a challenging 
task. A3395 is a multiple system formed of a t least three iden- 
tified clusters (see lTittlev & Henriksenl (l200ll) ). For the purpose 
of this paper we consider the two clusters most important in X- 
rays (SW and E), which can be clearly identified in Figure [3] 
The quality of the X-ray data is significantly poorer than in the 
case of A399-A401 (less than half the amount of exposure time) 
making it difficult to distinguish between background and dif- 
fuse X-rays in the intercluster region. Moreover, the clusters that 
form the system A3395 are clearly elongated (this is even clearer 
when studying XM M-Newton and Chan dra images of the clus- 
ter, see for example [Tittlev & Henr iksen (2001)) which adds an- 
other complication to the model. For the SZ part of the data, 
even though there seems to be a signal in the intercluster re- 
gion, the modelling of the SZ signal from A3395 is affected 
by the fact that Planck can not resolve the individual subsys- 
tems in A3395. Despite all these complications we are able to 
perform a basic fit to the SZ and X-ray data. The best-fit pa- 
rameters for the three pressure profile models are presented on 
the last three rows of tables [2] [3] and [4] respectively. We show 
in Figure |7] a) the tSZ Planck map for the A3391-A3395 pair, 
b) the best-fit GNFW2 model and c) the residuals after subtrac- 
tion of this model. We clearly observe in the residuals a sig- 
nificant excess of signal in the region near the A3395 system. 
Indeed, the integrated Compton parameter in the intercluster re- 
gion is Y = (4.5 + 0.7) 10~ 3 (arcminutes) 2 -of the same order of 
magnitude as that for the A399-A401 pair (see below). However, 
this residual is suspiciously close t o the NW component of the 
A3395 cluster that was discussed in iTittlev & Henriksenl d2001l) 
(and not included in our model) indicating that a far more com- 



Extra cluster 

In this case we assume the presence of a background cluster in 
the intercluster region of the pair. To simplify the modelling we 
consider that the cluster properties are fully defined by its red- 
shift and M500 mass which are the free parameters of the model. 
From M500 we compute R500 and use scaling relations to com- 
pute the temperature. We finally assume the pressure profile of 
the cluster is well descr ibed by a Un i versal profile and use the 
parameters obtained by lArnaud et al.l (T2010) for cool cores dis- 
cussed in the previous section. 



Filament-like structure 

In this case we assume that the intercluster region can be de- 
scribed by a tube-like filament oriented perpendicularly to the 
line of sight. We consider an isothermal filament of temperature 
T. For the electron density we consider two cases, one in which 
the density remains constant inside the tube and another in which 
the density falls from a maximum to zero in the border of the 
tube. The main physical parameters are the normalization den- 
sity, the temperature of the filament and its radius. For the latter 
the electron density profile is defined in the radial direction in 
cylindrical coordinates as follows 



n e {r) = 



n e (Q) 



(l+(t) 2 ) 5 



(7) 



with /3 = 4/3 and r c = 10'. The free parameters of the model are 
the temperature, T, and the central electron density, n e (Q). 
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Fig. 7. a) tSZ Compton parameter map, b) best-fit model of the clusters and c) residuals after subtraction of the best fit model for 
the A3391 and A3395 pair of clusters. 



6.2. The A399-A401 case 



From the three models considered in the previous section, we 
compute the residual signal both in tSZ and X-ray. The amount 
of signal in the intercluster region is then used to constrain the 
parameters of the extra cluster or the filament. In the case of an 
extra cluster behind the intercluster region we obtain the results 
in Figure [8] a) and c) where we trace the 1-D and (Figure 8b) 2- 
D likelihood contours at 68%, 95.5% and 99% confidence level 
for z and M500. Only clusters at high redshift z = 1 .95 and with 
a large mass M500 = 2.4 x 10 15 M o will be capable of produc- 
ing enough SZ signal with the strength of our Planck residual 
but with X-rays that do not exceed the ROSAT signal. These 
results are not consistent with existing constraints on structure 
formation that predicts that massive cl usters are formed at low 
redshifts (see for example Figure 2 in dHarrison & Colesll20ll . 



The best-fit parameters in the filament case are presented in 
Figure [9] where we show the 1-D and 2-D likelihood contours at 
68%, 95.5% and 99% confidence level for T and n e (0). We find 
that the temperature of the filament is (7.08 + 0.85) keV and the 
central electron density (3.72 ± 0.17) 10~ 4 cm" 3 . These results 
are co nsistent with previous estimates by ISakelliou & Ponrrianl 
(2004). It is interesting to look at the residuals for the best- 
fit clusters (we use here the GNFW2 model) and the filament 
model. Figure [TO] shows the 1-D longitudinal profiles and resid- 
uals for the tSZ effect and X-ray emission after subtraction of the 
clusters and filament models. We observe that the excess of tSZ 
effect in the intercluster region is well described by the filament 
model. The integrated Compton parameter in the intercluster re- 
gion is Y = (6.1 ± 0.7) 10~ 3 arcmin 2 . 
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Fig. 8. Constraints on the redshift and M500 for an extra cluster 
in the intercluster region of the A399-A401 system. 



7. Comparison with hydrodynamical simulations 

We have applied the full analysis described in the previous 
sections to hydrodynam ical simulations of a super-cluster-like 
region (iDolag et al.l l2006) that mimic our cluster pairs. Unlike 
the original simulations presented in IDolag etail (2006) 
this simulation was carried out including the treatment of 
radiative cooling, heating by a uniform UV background and 
star-formation feedback processes based on a sub-resolution 
model for the multiphase st ructure of the interstellar medium 
(Spri neel & Hernquisj 12003). The simulation also follows the 
pattern of metal production from the p ast history of cosmic star 
formation (Tornator e et al.l20 04. 2007). This is done by comput- 
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Fig. 10. a) tSZ longitudinal profile and b) X-ray longitudinal profile. Data Planck and ROSAT (black points), total model (red line), 
PXCC model only (light blue line) and filament model only (dark blue line). 
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Fig. 11. a) tSZ Compton parameter map (yx 10 6 ) for the Halo b in the hydrodynamical simulations, b) best-fit model for the clusters 
in the pair and c) residuals after subtraction of the latter. 
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Fig. 9. Constraints on the temperature and density of the filament 
in the intercluster region of the A399-A401 system. Note how 
a high temperature ~ 7 keV is also favored by the XMM data 
(SakelliojL&Ponman 20Ql 



ing the contributions from both Type-II and Type-la supernovae. 
Energy is fed back and metals are released gradually in time, 
according to the appropriate lifetimes of the different stellar 
populations. This treatment also includes, in a self-consistent 
way, the dependence of the gas cooling on the local metallicity. 
The fe edback scheme assumes a Salpeter initial mass function 
(IMF) (Salpeter 1955) and the parameters have been fixed to get 
a wind velocity of ss 480 km s ( - 1). 

We concentrated our analysis on a system of two merging 
clusters with characteristics similar to the A399-A401 system. 
In the simulation, a cluster with a virial mass of 6.5 x 10 14 M Q /h 
(Halo d) is merging with a cluster with a mass of 1 . 1 x 10 15 M Q /h 
(Halo b). At z = 0.07 the two systems are physically separated by 
4Mpc/h (2.9Mpc/h in the projection shown in figure 12) which 
compares remarkably well with the A399-A401 system, where 
the projected separation of the systems is ~ 3 Mpc. More de- 
tails on the appearan ce and the geomet ry of the simulated sys- 
tem can be found in iDolag et aD (120061) . Figure QT| shows the 
Compton parameter maps (multiplied by 10 6 ). After subtracting 
the best-fit model for the clusters (Figure 1 lb) we fit the resid- 
uals (Figure 11c)) in the intercluster region with the filament 
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model described above. For the best-fit model of the filament the 
integrated Compton parameter is Y — (3.5 + 0.7) x 1CT 3 arcmin 2 
and is in good agreement with the input Yj nput = 3.2 x 10~ 3 
arcmin 2 . It is also important to notice that the error bars obtained 
are consistent with those obtained for the cluster pair A399- 
A401. 

We trace the origin of the particles causing the signal within 
the region between the two clusters. We therefore distinguish 
between particles belonging to the main part of one of the two 
galaxy clusters (e.g. within 0.5 R V j> 13) and the particles belong- 
ing to the filament, selected as those particles within a cylinder 
of IMpc/h in diameter between the two clusters and lying out- 
side 0.5 R V! - r . Additionally, we also define particles in the central 
part of the filament as those within a cylinder of 0.5 Mpc/h in 
diameter between the two clusters and lying outside 0.75 R V(r . 
Figure [12] a-c shows the different location of the particles color 
coded in blue and green for particles belonging to the clusters, 
red for particles in the filament, and cyan for those in the inner 
part of the filament. The yellow circles mark the virial radius of 
the two clusters. Going back in time to z=0.1, z=0.2 and z=0.3 
(corresponding to 0.4, 1.5 and 2.5 Gyr) reveals the spatial origin 
of these particles. We can observe that according to the origin 
of these particles, the intercluser region at z=0.07 is populated 
by two separate components. One component is coming from 
material that belonged to the outer atmosphere of one of the two 
clusters at z=0.3, and the other one is coming from a very diffuse 
structure, outside but connecting both clusters. This structure has 
a sheet-like geometry at early times. When we trace back the ori- 
gin of the particles in the very central part of the filament, most 
the particles were further than 0.5 R V! > from one of the main 
clusters within the last 2.5 Gyr. This fact is made more evident 
in Figure[13] which shows the spatial distribution of the particles 
of the different components more quantitatively. Here we sorted 
each particle by its minimum distance to either cluster. When we 
concentrate our attention on the filament's core (light blue par- 
ticles), only a very small fraction of the particles in this region 
were closer than 0.5 R v ,„- to one of the two main halos at z=0.3. 
More than half of them are coming from regions with 0.5 R v „ < 
R < R l7 > and a non-negligible fraction is coming from regions 
beyond R Vi >. 

This simulation shows that the intercluster region might be 
populated by a mix of components from both the clusters and the 
intergalactic medium. 

8. Discussion 

The presence of significant SZ signal between pairs of clusters 
combined with the lack of significant X-ray signal suggestes 
that merging events may have occurred, leaving filaments of 
matter between these interacting pairs. Merging of the clusters 
would contribute to the increased pressure that boosts the 
SZ signal. These results are in qualitative agreement with 
hydrodynamical simulations in which the intercluster region 
betwee n an interacting p air of clusters can show significant SZ 
signal dDolag et alj|20 06). 

Regarding the cluster pair A399-A401, earlier studies 
showed that the intercluster medium in this pair is compresse d 
due to the mergin g process (Fujita e t al. | Fuiita et al.l fl996), 
Sakelliou et al. ISakelliou & Ponmanl (2004)). The increased 
pressure would increase the SZ signal. (AJcahori & Yosh ikawal 

1 We define R n > based on the over-density from the top hat spherical 
collapse model. 



(120081) 1 study the non-equilibrium state of this system with sim- 
ulated data and find that there might exist significant shock lay- 
ers at the edge of the linked region between the clusters that 
could explain a boost in the SZ signal. Also, Suzaku observa- 
tions found that the inte rcluster medium h as a relatively high 
metallicity of 0.2 solar (Fuj ita et al.l (120081) ). These works es- 
timated that the filame ntary bridge would have an electron den- 
sity of n e ~ 10~ 4 ctr T 3 (Fujita et al.l dT996l) . ISakelliou & Ponmanl 
(2004. lFtffita1st"aLl (120081) 1. 

Given that the angular separation between the two clusters is 
~ 3 Mpc, and that both clusters are at slightly different redshifts 
(0.0724 and 0.0737), the true separation between the clusters is 
~ 5 Mpc. That is, the clusters are separated by more than their 
respective virial radii. Consequently we can conclude that the 
signal (or at least a significant percentage of it) seen by Planck in 
the intercluster region corresponds to baryons outside the clus- 
ters. 

Our results show that there is evidence for a filamentary 
structure between the pair A399-A401 and outside the clusters 
but also opens some questions about the origin of this gas. The 
uncertainties in the estimation of the cluster contributions around 
their virial radii makes it difficult to discriminate between the 
different scenarios. In a pre-merger scenario a filament could 
have been trapped in between the clusters, with its material being 
reprocessed and compressed as the clusters approach each other. 
In a post-merger scenario (meaning by post-merger not a direct 
crossing of the cluster cores but a gravitational interaction as the 
two approaching clusters orbit their common centre of mass), 
the intercluster signal could be just the result of the overlapping 
tails of the disturbed clusters. 

These clusters could be elongated due to their gravitational 
interaction or a bridge of matter has been formed between the 
clusters after an interaction. The fits to SZ and X-Ray data 
reveals that there are some apparently conflicting results, that 
make it harder to reconcile the SZ and X-ray data with spherical, 
standard models. A good e xample is the best fitting model of 
ISakelliou & Ponmanl d2004l) that over-predicts the SZ data. This 
could be, for instance, an indication that spherical models might 
not be the most appropriate for describing these clusters. 

Non-sphericity is expected (and indeed observed) in X-rays 
and could introduce corrections to the best fitting models, partic- 
ularly if the elongations of the clusters point towards each other 
(as it is the case for A401 from the X-ray images). In this case, 
the gas density in the intercluster regio n would be enhanced due 
to this elongation. ISereno et al.l ((2006) shows how clusters seem 
to favor prolate geometries. A prolate model would increase the 
X-ray signal for the same number of electrons (or SZ), or vice 
versa, the same X-ray signal requires fewer electrons (or SZ). 
We have estimated that a ratio of ~ 1.3 between the axis might 
be sufficient to reduce the SZ signal by ~ 30% with a fixed 
X-ray observation (i.e, fixing the other two axis and the total 
X-ray signal). 

Clumpiness is another factor that might be introducing a 
bias in the models tha t best subtract the cluster components. 
Mathiese n et al.l (1 19991) use simulations to estimate the mean 
mass-weighted clumping factor C -< n 2 > / < n > 2 . They 
found typical values for C between 1.3 and 1.4 within a density 
contrast of 500. Since for an isothermal model C can be seen as 
the ratio of X-ray to SZ signal (squared), as in the case of prolate 
clusters discussed above, clumpiness of the gas would increase 
the X-ray signal for the same number of electrons (or SZ) or 
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Fig. 12. Spatial projections (along x,y and z axis) of the particles in the simulation. The rows show the simulation at z=0.07, z=0.1, 
z=0.2 and z=0.3 respectively. We colour code the particles depending on their location: clusters (green/blue), the filamentary region 
in between (red) and the central part of this filament (cyan). See text for a more detailed definition. Additionally the yellow circles 
mark the virial radius of the two clusters. 



vice-versa, the same X-ray signal would require fewer electrons 
(or SZ). 

One cluster will have an effect on the other cluster and in the 
intercluster region through its gravitational potential. Figure FBI 
shows the resulting gravitational field across the line intersect- 
ing the two cluster centers. The gravitational field is computed 



from the hydrostatic equilibrium equation assuming a profile 
of the gas density (and a constant temperature). We use one 
of our best fitting models described earlier (universal profile) 
for this purpose. The intercluster region shows a significant 
enhancement in the gravitational potential. The effect is more 
significant in A399 where the level in the intercluster region is 
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Fig. 13. Evolution of the radial distribution of the particles with respect to the centre of the closest of the two clusters. Colour code 
is the same as in figure[T2] The dashed and dot-dashed lines correspond to the blue and green particles respectively (clusters) while 
the dotted and solid lines orrespond to the red and cyan particles (filament). 



higher with respect to the maximum acceleration at the centre 
of the cluster. Therefore we should expect the gravitational 
attraction experienced by the gas in the clusters to be stronger 
(compared with the peak of the potential) over A399 than over 
A401. The gravitational field in the intercluster region increases 
as the two clusters approach creating a pulling effect over the 
gas of the clusters towards this region. In this scenario, since 
the gas would be moving from gravitational fields with similar 
intensity it would not undergo an adiabatic expansion and hence 
would retain its original temperature. This would explain the 
high temperature found in the intercluster region. This scenario 
would also explain the high metallicity in the intercluster 
area (Z=0.2, Fujita el al. 2007, although the constraints in the 
metallicity are not very strong and should be take with caution). 
If the gas in the intercluster region was originally in a filament 
it would be difficult to explain such a high metallicity (if 
confirmed) but not if the gas comes originally from the clusters. 

The gravitational pull of the intercluster region over the outer 
parts of the cluster may possibily explain how the intercluster 
region can be partially populated with metal-rich gas from the 
clusters. As the simulation presented in section [7] shows, this is 
a reasonable scenario however, it also shows that the interclus- 
ter region can be populated by intergalactic material. Thus the 
intercluster region might be a mix of cluster and intergalactic 
material. 

9. Conclusions. 

Using the Planck and ROSAT data, we have studied the tSZ 
and X-ray maps of 25 pairs of clusters of galaxies. After mod- 
elling (assuming a spherical symmetric model) and subtracting 
the contribution of each individual cluster, we have detected sig- 
nificant tSZ residuals in at least two of these pairs : A399-A401 




2 4 

Mpc 



Fig. 14. Gravitational acceleration across the line connecting the 
two cluster centers (A401 left and A399 right. The accelera- 
tion is computed assuming hydrostatic equilibrium. The model 
GNFW1 was assumed for the computation of the hydrostatic 
equilibrium equation 

and A3391-A3395. In the case of the A399-A401 pair, these 
residuals are compatible with an intercluster filamen t of hot, 7 
keV (in agreement with Sakellio u & Ponma 3 d2Q04)), and dif- 
fuse, 3.7 x 10 -4 cm -3 , gas connecting the two clusters. A chance 
coincidence of a background cluster is ruled out for canonical 
scaling relations as it would have to be a very massive cluster 
(M 500 = 2.4 x 10 15 M Q ) at a very high redshift (z = 1 .9). The sig- 
nal detected by Planck is significant independently of the cluster 
model. Hydrodynamical simulations show that the intercluster 
signal in A399-A401 is compatible with an scenario where the 
intercluster region is populated with a mixture of material from 
the clusters and the intergalactic medium indicating that there 
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might be a bridge of matter connecting the two clusters. If the 
measured signal of the merging cluster pair can be interpreted in 
terms of spherically symmetric individual clusters, evidence re- 
mains for an intercluster SZ signal detected by Planck. It is con- 
sistent with simulated data and may constitute the first detection 
of the tSZ effect between clusters. Under this interpretation, the 
signal is unambiguous in the sense that it is detected with high 
significance (as shown by Figure 2), is not due to any known arti- 
fact, is well resolved by Planck, and is located in what one would 
consider the external region of a standard cluster. The exact in- 
terpretation of the origin of the signal is more open to specula- 
tion. The analysis presented in this work shows the potential of 
Planck data to study these yet unexplored regions. Better angular 
resolution observations of the tSZ would improve the modelling 
of the clusters and reduce the uncertainties in the estimation of 
the signal excess in the intercluster region. 
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